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ABSTRACT 

We obtain a bilocal classical field theory as the large limit of the chiral Gross- 
Neveu (or non-abelian Thirring) model. Exact classical solutions that describe topological 
solitons are obtained. It is shown that their mass spectrum agrees with the large N limit 
of the spectrum of the chiral Gross-Neveu model. 
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We will study the large N limit of the two- dimensional fermionic system defined by 
the Lagrangian 

L = q'[-i^,d^]q^ + ^q^^^q^q^^^q,. (1) 

Here i = 1, - ■ ■ N. This model is a non-abelian version of the Thirring model. It is also 
known as the chiral Gross-Neveu [1] model since the interaction can be rewritten in the 
following form using a Fierz identity. 

L = q'[-il^d^']qi - ^Wqiq^qj - q'l^m'l^qj]- (2) 

The precise definition of the model requires a renormalization of the coupling constant. 
This theory is asymptotically free, the dimensionless coupling constant is replaced by 
a dimensional constant as the true parameter upon renormalization. This model has been 
solved exactly [2], [3] by Bethe ansatz methods; i.e., its spectrum and the S'-matrix are 
known explicitly. In the limit as ^ oo, it should tend to an exactly integrable classical 
theory. This classical theory will require a renormalization in order to be well-defined ( the 
beta-function does not vanish in the large N limit). Still, we should be able to understand 
the large limit of the mass spectrum in terms of the classical solutions. Such a direct 
understanding of the large N limit of this type of theories has not been achieved yet, as 
pointed out in [4]. Berezin [5] has given a bi-local formulation of the large limit of 
our type of model; however, the mass spectrum was not obtained in that reference. What 
amounts to the linear approximation to our theory has been studied in a recent paper [6] . 

In this paper we will obtain the classical theory corresponding to the large N limit 
of this model as well as obtain its mass spectrum, in agreement with the known exact 
solutions. We will apply the bilocal bosonization method [7], [8], [9] as well as the 
non-perturbative renormalization methods [10] developed in previous papers to this 
model. Having established the validity of these methods in this exactly solvable context, 
these methods can in the future be used in more realistic models such as spherical QCD 
[8]. In fact, one of the lessons of this paper is how to define classical field theories which 
require a renormalization. We believe this is the first example of exact soliton solutions in 
a bilocal field theory. We also discover some new phenomena such as a topological charge 
with a continuous range of values. 

The large N limit of our model can be studied by the technique of summing an 
appropriate class of Feynman diagrams. As is clear from other examples, [11] this 
will produce a free field theory of small oscillations around the vacuum, the analogues of 
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mesons in two dimensional QCD. However, it is clear from those examples that the true 
large N limit is a highly non-linear classical theory. In particular it can have classical 
solutions which are topological solitons and hence are large deviations from the vacuum. 
The bi-local bosonization methods of Refs. [7] can construct the complete large N limit 
and describe even these topological solitons [12]. ( In 2dqcd they are the baryons.) 
Indeed we will see that in the case of the chiral Gross-Neveu model, all the particles in the 
spectrum are topological solitons. A surprise is that the topologically conserved quantum 
number can take continuous values between and 27r. A completely precise mathematical 
description of such a topological invariant, with a continuous range, requires the theory of 
projectors on von Neumann algebras and is beyond the scope of this paper. The technical 
apparatus necessary for this seems to exist already [13]. 

Let us begin by understanding the symmetries of the above lagrangian. First of all 
it has a global symmetry under U{N) x U{N), the first (second) U{N) acting on the left 
(right) components of q. This non-abelian part of this chiral symmetry is unbroken and 
the particles of theory transform under the completely anti- symmetric tensor representa- 
tions of SU{N) X SU{N). ( This is unlike in 2dqcd, where all particles are in the trivial 
representation of SU (N) due to confinement.) If r = 1, ■ ■ ■ , — 1 is the rank of the tensor, 
the particles transform under the representation (r, r) of SU (N) x SU (N) and have masses 
[3], [2] 

mr = mi- — —. (3) 
sm^ 

The quantum number r is also a conserved quantity of the fermionic lagrangian, corre- 
sponding to the symmetry group Z^- We might understand this discrete symmetry by 
considering the U{1) sub-group of the global symmetry corresponding to 

Qi ^ e'^Qi. (4) 

If a is an integer multiple of this transformation is in fact a particular element in 
SU{N) x SU{N). This Zn subgroup may viewed as measuring the fractional part of 
fermion number, if is assigned a fermion number of ( This is analogous to the 'triality' 
of SU{3) representations.) The part of fermion number which is integral can be viewed as 
a separate conserved quantity, whose U{1) symmetry group has no overlap with SU{N) x 
SU{N). We will thus decompose fermion number into two conserved quantities, one taking 
integer values on each completely occupied one-particle state (which we will call baryon 
number B in analogy with 2dqcd) and another, 9 taking values 0, 27r;^, ^tt-^, • • • , 27r[l — jj^] 
corresponding to Zn. 
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In the Dirac sea, any one-particle state that is completely filled will contain N fermions 
and will make no contribution to the conserved charge 9 of this symmetry. If a state is 
partially filled, its contribution will be equal to the ratio of the number of occupied states to 
the number of available states, which is an integer multiple of jj^. The conserved quantity 
^ is precisely this fractional part of the baryon number. Obviously, the 'fractional part' 
of the fermion number can add up to give an integer if we sum over many states. Less 
obvious is that the integral part of the fermion number can add up to a fractional value 
when summed over an infinite number of states; this has to do with regularizations that 
are necessary to make such a sum well-defined [14] . 

There are several possible ways of taking the limit N ^ oo. For example, we could 
take the limit keeping the mass mi of the lightest particle fixed. Then, the spectrum would 
consist of two infinite towers of particles of equally spaced masses; one tower corresponding 
to r = 1, ■ ■ ■ , ^ and the other tor = A^ — 1,A^ — 2,---^ + l. ( This is for N even; there 
is a similar expression for odd.) A more interesting large limit will be obtained if we 
take the limit keeping the mass /j, of the heaviest particle fixed. Then the mass differences 
between particles tends to zero in the large limit and the spectrum merges into a 
continuum. (Such limits of the exact solutions have been studied before in a somewhat 
different context. [15].) Also, in this limit, the variable 6 takes a continuous range of 
values between and 27r so that the mass spectrum becomes the continuous set, 

mo — /Ltsin -, for < 6* < 27r. (5) 

We will try to understand this limit as a classical field theory. We will see that Q corresponds 
to a topologically conserved quantity with a continuous range of values between and 27r 
in the classical theory. 

First we will give a heuristic derivation of the classical theory corresponding to the 
large N limit of the Thirring model. Then we will give a more rigorous ah initio definition 
of this classical theory and find some static solutions to it. We regard the following 
discussion of the large N limit only as motivation for the study of this classical theory. A 
more rigorous justification would be in terms of the quantization of the classical theory, 
along lines discussed in [7] . will play the role of Ti in this quantization. 

We will study the theory in hamiltonian form, in which the field operators satisfy the 
equal time anti-commutation relations 

[^l* (^), 9,73(y)]+ = ^c.^^]K^ - y)- (6) 
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a, /3 = 1, 2 are spin indices, which we will find convenient to suppress usually. If we ignore 
delicate issues of renormalization, the hamiltonian can be expressed in terms of the color 
singlet bilinear M{x, y) = —jj^ '■ (li{x)q^'^{y) : ( here, M{x, y) is a 2 x 2 matrix in spin space 
and : : denotes normal ordering with respect to the free hamiltonian): 



H 
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tr(-i75)[^^^^^],^^ - |-tr[M^(x,x) -75M(x,x)75M(x,x)] 



dx. (7) 



Conventional large N arguments can be used to show that the observables M{x, y) will 
have quantum fluctuations of order so that their time evolution will be described by 
classical equations of motion. (One way to see this is to note that their commutators are 
of order [7] . They form a representation of the infinite dimensional unitary Lie algebra, 
also called W^q algebra in the context of matrix models.) In a theory such as 2dqcd which 
has color confinement, these would be a complete set of observables. However, the physical 
states of our theory transform under non-trivial representations of color, so that these are 
not a complete set of observables. 

We can now see that M{x,y) along with the charges of the SU{N) x SU{N) sym- 
metry are in fact a complete set of observables. It is best to understand this fact in a 
regularized context in which the position variables are allowed to take only a finite number 
( say K) of values. The complete set of bilinears ^i^{x,y) = — i[g*"(x), gti/3(y)] form the 
Lie algebra U_{2KN) under commutation. The Fermionic Fock space carries an irreducible 
representation of this Lie algebra. This Lie algebra contains SU (N) x SU{N) x U_{K) as 
a subalgebra. The generators of this subalgebra are the charges of the global symmetry 
SU (N) X SU (N) and the operators M{x, y) themselves. The key fact is that the repre- 
sentation of U {2KN) on the fermionic Fock space remains irreducible with respect to this 
subalgebra. ( We omit a detailed proof in the interest of brevity.) Thus any operator 
that commutes with the charges of the SU{N) x SU{N) symmetry as well as M{x,y) is 
a multiple of the identity: these together form a complete set of observables. Thus the 
dynamics of the theory reduces to a classical dynamics for M{x,y)] the generators of the 
global symmetry never become classical, but they have trivial dynamics, being conserved 
quantities. Thus we will for the most part concentrate on the variables M{x,y). These 
arguments can probably be made rigorous in the case were x,y have infinite range, but 
the techniques required probably involve use of C* or von Neumann algebras. We will not 
attempt it in this paper. 

Any classical theory is described by a phase space ( manifold of allowed initial con- 
ditions), a symplectic structure on this manifold ( or, Poisson brackets of a complete set 



5 



of dynamical variables) and a hamiltonian. In some cases, the phase space is described as 
the set of solutions to some set of constraints satisfied by some dynamical variable, rather 
than directly in terms of a co-ordinate system. 

In our classical theory, the dynamical variables are self- adjoint operators M on the 
complex Hilbert space L'^{R, C^). ( This can be thought of the one-particle Hilbert space 
of a Dirac fermion in 1-|-1 dimensional space- -time.) M can be described in terms of its 
kernel, 

Mu{x) = j M{x,y)u{y)dy (8) 

where M{x,y) is a 2 x 2 matrix valued function ( in general distribution) of the two 
variables x,y. An equivalent description is in terms of the 'symbol', M{x,p), which is a 
Fourier transform of the kernel with respect to the relative co-ordinate: 

M{x,p) = J M{x+^,x-^)e-'Pydy. (9) 

Since M is self-adjoint, M{x,p) is real valued. Spatial translation acts as follows on these 
variables: M{x,y) M{x + a,y + a) and M{x,p) — > M{x + a,p). We can regard M{x,p) 
as an infinite component classical field, labelled by a continuous internal index p. Such 
bilocal field variables have been useful in formulating the large Nc limit of models for QCD. 

In spite of the fact that operators on Hilbert spaces make an appearance, we emphasize 
that ours is still a classical theory, the dynamical variable of which just happens to be an 
infinite dimensional 'matrix' M{x,y). The Poisson brackets satisfied by these variables 
are, 

{M^{x,y),M^{z,u)} = d^d{y - z)[Mi{x,u) + eUx,u)] 

-Si5{u-x)[M^iz,y) + e^^^{z,y)]. 
Here a, /3 are spin indices which we usually suppress. Also, 

eLi^, y) = iL j ^e^^(--^) sgn (p) (10) 

is the kernel of an operator eo whose square is one. It has the physical meaning of being 
sign of the massless Dirac operator; the eigenspace of eo with eigenvalue 1 (or -1) is the 
space of positive (or non-positive ) energy one-particle states. 

This algebra is the central extension of the unitary Lie algebra, called the Lundberg- 
Kac-Petersen [16] extension. This algebra has also acquired the name Wqo algebra 
in the context of matrix models. The above Poisson brackets are just the commutation 
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relations of the operators M{x,y) introduced previously, except that a factor of has 
been removed. This is the appropriate prescription as plays the role of h in our classical 
limit. The central term proportional to eo appears because of the normal ordering of M 
in the quantum theory. 

We will impose some conditions on the asymptotic behaviour of M{x,p) as \p\ oo. 
These follow from the asymptotic behaviour of the matrix elements of 

in free fermion theory. Due to asymptotic freedom, this is the same as the behaviour of 
the interacting theory. However in the classical theory they are to be viewed as postulates. 

We will impose 

Maix,p)^Oij^) Mod{x,p)r^Oi^) (11) 

\p\ \P\ 

where Md{x,p) (or, Modix,p)) is the diagonal (or, off-diagonal) part of the 2x2 matrix 
M{x,p) in a basis where 75 is diagonal. It will not be possible to absorb the central term 
into the definition of M{x, y) without violating these conditions; this is what makes the 
central extension non-trivial. The above commutation relations then define a topological 
Lie algebra, the topology being defined by the norm implicit in this asymptotic behaviour: 

||M|| = sup[(/ + l)\Ma{x,p)\] + sup[(|p| + l)\Mod{x,p)\]. (12) 

x,p x,p 

The phase space of the theory is a co-adjoint orbit of this Lie algebra ( or rather the 
corresponding Lie group, which is a dense subgroup of the restricted Unitary group Ui{H) 
of Segal [17]; we refer to [7] for a more detailed description). In the case of 2dqcd, this 
orbit was a Grassmannian, which is defined by a quadratic equation satisfied by M: 

[M + eo]' = l. (13) 

The co-adjoint orbits of the restricted unitary group are known; they are infinite 
dimensional analogues of the familiar fiag manifolds. ( This theory is not necessary in the 
following, except as motivation for some definitions.) We will require that M satisfy the 
above constraint except for a finite dimensional block. More precisely,that 

[M + eof -1 is finite rank. (14) 
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For the static solutions we are interested in, this block is in fact one-dimensional. The 
meaning of this condition is that we allow for fermion states to be completely filled or 
completely unfilled except for a finite dimensional block of states which may be only 
partially filled. If in fact (M + eo)^ = 1 all the states in the theory would be either 
completely filled or completely unfilled and therefore singlets under the global SU{N) 
symmetry. [18] This is too strong a condition in our case. In fact the constraint that 
(M + eo)^ — 1 be finite rank will enforce that 

M\x,p)^ -sgn{p)['y5,M{x,p)]+ (15) 

or that 

M^{x,p) + Mld{x,p) ~ - sgn {p)[^^,Mod{x,p)] + . (16) 
These asymptotic conditions will be useful later. 

The hamiltonian of the theory can be obtained by rewriting the regularized fermionic 
hamiltonian in terms of the bilinears. It is convenient to write it in terms of M[x,p), for 
then a cut-off in the range of p provides a natural regularization. After some manipulations 
we can bring the hamiltonian to the form, 

Ea{M) = fdx[[ ^tr[75pM(x,p) + lv{x)M{x,p)] - C]. (17) 
J J\v\<A 27r 2 



Here, 



V{x) = 2g\K)l MUx,p)^ (18) 

and 5'^(A) is a coupling constant which will be picked such that the theory is well-defined 
( see below) . Also C is picked such that the vacuum solution has zero energy density. The 
equations of motion are anyway independent of C. 



It is clear from the asymptotic behaviour of Mod{x,p) that , if g'^{A) V{x) 



log A ' 

will be independent of A asymptotically. In fact, this will turn out to be a self- consistent 
choice for 5'^(A). If moreover, (/f = f , the quantity 



tr[75pMd(a;,p) + ^V{x)Mod{x,p)] (19) 

vanishes faster than -j^ . ( Each term in the trace goes like but the leading contributions 
cancel.) This means that the energy can be defined by taking the limit A — > oo, 

E{M) = j dx[j ^tr[75pM(a;,j9) + ^V{x)M{x,p)] - C] ; (20) 
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a regularization is only necessary in the above expression for V{x) in terms of M{x,p). 



The equations of motion that follow from this hamiltonian and the above commuta- 
tions relations are, ^^g^^'^'*'* = {E{M), M}. The Poisson bracket relations being those of 
a imitary Lie algebra, the r.h.s. can be written in terms of the commutator of operators. 
If we use operator notation, we get, 

^^ = [H{M),M + eo] (21) 
where the operator H is defined to be the derivative of E with respect to M, 

nm = ^. (22) 

Explicitly, H{M) is a differential operator, 

H{M) = -i^,-^ + Vix) (23) 

where V{x) is related to M by the equation given earlier. Static solutions must therefore 
satisfy the nonlinear equation 

[M + eo,H{M)] = 0. (24) 

Since we are studying the large N limit of a theory which is exactly solvable for every 
A'', there must be a way to solve these equations as well: the above classical dynamical 
system must be integrable. We do not attempt to demonstrate the exact integrability of 
this system here; instead we will obtain a family of classical static solutions that describe 
a kind of topological soliton. We will show that the mass spectrum of this soliton agrees 
with the known large N limit of the mass spectrum of the Thirring model. 

A simple solution to the above equations of motion is, 

M{x,p) + eo{p) = ^^^2+^2) (25) 

with V{x) = mjo- The relation between V{x) and M is satisfied if 

2g\A) [ —--^ ^ = 1. (26) 

This solution describes the vacuum, being translation invariant. From now on we will 
assume that g'^(A) is given by the above equation. We have traded the dimensionless 
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coupling constant g for the dimensional constant m. The constant C is now fixed so that 
this solution has zero energy: 



(27) 



which is a convergent integral. 



In fact the vacuum is degenerate; we could have replaced it by a chirally rotated 
solution, 

V(x) = e-*T'^tm7oe'^^^ Mix.p) + €o(p) = e'^^^i ^ff^^ e'^^'^ (28) 

which would have the same energy ( zero). This suggests the possibility of more general 
solutions which depend on x and tend to two different vacua as |a;| — > ±00. Such solutions 
would describe topological solitons of the bi-local theory. Notice that there is a continuous 
infinity of vacua and thus the solitons are parametrized by an angle ^, measuring the 
difference between the directions of the vacua at infinity. We do not know of a systematic 
way to search for such static solutions of the above classical equations for M. Instead 
we will produce a solution by a guess inspired by the theory of solitons of the non-linear 
Schrodinger equation. 

Now, H is the Dirac operator of a massive spin \ particle coupled to an external 
scalar field Vix). H has both a point ( discrete) spectrum corresponding to bound states 
and a continuous spectrum corresponding to scattering states. The discrete spectrum will 
have eigenvalues in the range —m < X < m and will have eigenfunctions which are square 
integrable.The continuous spectrum ( 'scattering states') have eigenvalues with A < — m 
or A > m. There is no normalizable eigenvector corresponding to these; however, it is still 
meaningful to speak of the projection operator Pj to the subspace with eigenvalues within 
some interval / of the real axis. This spectral projection operator Pj can be written as 

P, = / ^P(A) (29) 

where the 'spectral density' p(A) is in general some operator- valued distribution. The 
contribution of a bound state will involve a delta-function in A while that of the contin- 
uum will be some continuous function in A. In terms of the spectral density, we have a 
decomposition, 

/r/A 
-Ap(A) (30) 
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which is the generahzation of the famihar decomposition of a matrix into eigenvalues and 
eigenvectors. We can determine p(A) as the discontinuity the resolvent -R(A) = {H — X)~^ 
of H across the real axis: 

p(A) = lim ^ [R{X + ie) - R{X - ie)] . (31) 

The resolvent in turn can be determined in terms of the Jost solutions [19] of scattering 
theory. Thus if we have an explicit expression for V{x), we can in principle solve the 
scattering problem for H and obtain p(A). 

The static equation of motion implies that H and M + eo commute; one way to satisfy 
this is for M + eo a-nd H to be simultaneously diagonal. Gcncrically there will be no 
degeneracies and this will the only solution. Thus we should expect there to be a real 
valued function a{X) such that 

M + eo = J a(A)p(A)£. (32) 

This function must approach ±1 as A — ±oo to satisfy the asymptotic condition. Its 
value in a region not contained in the spectrum of H will not matter to M, since p{X) will 
vanish there. Any discrete eigenvalue of H will be contained in an interval in which it is 
the only element of the spectrum. Thus we can choose to extend (j( A) to be a constant in 
a neighborhood of every discrete eigenvalue. This choice will turn out to be useful later. 

We can regard M{x, y) as describing the self-consistent way to fill the energy states 
of H with fermions.In fact, is the filling fraction, the number of states that are 

occupied as a fraction of the total available ( which is N). For large positive A, this must 
go to zero and for large negative A it must go to one; there can be a finite number of states 
with a fractional filling factor. 

Now in our problem, V{x) and M are related by the self-consistency condition. 

We have already fixed g'^{A) from the vacuum solution, so that the r.h.s. is logarithmically 
divergent. We will now show that for any V{x), the r.h.s. also contains a logarithmically 
divergent piece; after cancellation of this piece, we will get a convergent equation. 

First, note that by a change of variable from p to A, 

1 r dXQiX^ -m'^) 2 ^2^ / 
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Here denotes the step function and 

k{X) = V(A2 - m^). 
Now, the consistency condition can be written as. 



(35) 



dX 
2^ 



J\p\ 



<A 



dp 
2^ 







Now, for large A, the leading contribution of the second term will come from the region of 
large \p\ . In the limit of large p, the asymptotic behaviour of p( A; , p) will be determined 
by that of the resolvent symbol ^(A; x,p). Furthermore, the resolvent symbol for large \p\ 
is, 

R{X;x,p)r.[^,p + V{x)-X]-^ (36) 



since the WKB approximation applies in this case, for smooth V{x). Thus we find. 



where. 



and 



p(A; x,p) ~ 5{X - A+)n+ + 6{X - A_)n. 



A± = ±V(/ + I^^P) 



n± 





V 


2 


\v 


2 + (A 





1 



>^±-p 



A±-P {X±-pf 



(37) 
(38) 
(39) 



V \v\'' 

In the above, we have assumed without loss of generality that 



V{x) 



v{x) 
v*{x) 



(40) 



and we have sometimes omitted the x,p dependence of v, X±, Ilpm etc. for simplicity. With 
this asymptotic expression, it can be verified that the divergent terms in the consistency 
condition cancel out; so we can take the limit as A — > oo to get 



/ 



dX 
2^ 



a{X)pod{x,x; X)-0{X'-m^)^ 



= 0. 



(41) 



We have used here, / ^pod{x,p) = pod{x,x). Although each term separately would be 
divergent, the quantity in the square brackets has a finite integral over A. (This can also 
be verified using the explicit forms given below.) 

It remains now to solve the above nonlinear integral equation for V{x) and a-(A). We 
will make a guess and see if it in fact satisfies the equation. There must be a systematic 
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method using the inverse scattering theory, [19] but we will not attempt to develop this 
here. 



We propose the ansatz, 



TTl 

V{x) = mjo + -j-^ (70 - 70) (42) 
r]{x) 



where, 

a 

r]{x) = 1 + e-"'' , v = 2m sin - (43) 

and 

7, = Q-i(^)7oQ(^), Q(^) = e^75e/2_ (44) 
As X — > —00, it tends to the vacuum solution, V{x) — > m^Q. As a: ^ 00, it tends to 

—iga^d i^Q^ 

another vacuum solution, differing by a chiral rotation ■.V{x) — > e 2 771706^". Thus 
this ansatz describes a sort of topological soliton of our theory. 

This ansatz for V{x) is the well-known [19] reflectionless potential of the Dirac op- 
erator, which is known to be a soliton of the nonlinear Schrodinger equation. We will 
use these potentials to produce solitons of our bi-local field theory. There is no direct 
relationship between the nonlinear Schrodinger equation and our bi-local theory: indeed 
our theory is relativistically invariant while the nonlinear Schrodinger equation is invariant 
under Galilean transformations. Yet, the static solution of both systems involve the same 
reflectionless potential!. 

The point spectrum oi H = —^75^ + V{x) consists of one normalizable eigenstate 
('bound state'). 



with eigenvalue Xb = ?77,cos|. Note that — m < Xb < m, so that the bound state is in 
the 'gap', the set of values that are forbidden as eigenvalues of the free Dirac operator. 
The continuum eigenfunctions ( scattering solutions) can be written in terms of the Jost 
functions, [19] and we can obtain the answer for p(A) explicitly. We omit the computations 
and just display the answer *: 

p{x, y; A) = 27r5(A - Ab)V'b(x)V'b(2/)+ 



2k \(3{x)a*{y) + a*{x)l3{y) a*{x)a{y) + /3{x)l3*{y) 



from now on we use a basis in which 75 = ( ^ ^1 J > 7o ^ * 



-1 y ' \i 
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Here, 



: <t> + 



Moreover, is a sort of scattering phase shift and is determined by the transcendental 
equation 

sm— — = sin^^-— . (47) 



We only need the special case pod{x,x) to verify the equation of motion. We put in 
as part of the ansatz, cr(A) = sgn (A) in the continuum, |A| > m. After some calculations 
the consistency condition becomes. 



I 



"'^^/x2 2N'^™r/-, e""-!, ,^ e* 2 - 1 2 - 1 



Usina; -%r^ = -rr 3^ we see that the terms are proportional to either or -^f-^ 

( Terms independent of 77(0;) cancel out.) The terms proportional to give, after some 
simplifications, 

o-(AB)msin 2 = ~ / 2^"fc ^^'^ ~ m^)4sin sin (48) 

The terms proportional to happen to give the same equation. With (j(Ab) given by 
the above equation we have a solution to the static equations of motion. 

The expression for o"(Ab) can be simplified further. If we hold 6 fixed, the scattering 
angle <p can be thought of as a function of A: either using the previous transcendental 
equation, or the alternative forms, 

6-0 9 k + X 6 k + X + m 

cot — - — = cot - H 5-, tan— = - tan-. (49) 

4 2 msinf 4 k + X - m 4 

If we change the variable in the integral from A to (j), using the differential of the first one 
of the above equations 

d(j) A dX 

we will get 



(7(AB)msin^ = -4 / ^[1 + ^]G(A2 - m^) sin^ 
2 J 27r k 



4 



1 sin 
— m— — 
27r 2 



\=—m n\=oo 

d(j)+ dcp 

= — 00 J \=m 
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Since = —9, 0, 2tt, 6 respectively at A = — oo, —m, m, oo, we have 

a{XB) = (50) 

TT 

To summarize, a static solution of the classical equations of motion is determined by- 
two real functions cr(A) and V{x) satisfying the self-consistency equation 



2^ 



2 ^2.V{X) 



cr(A)pod(a^,a;; A) - ©(A - m ) 



(51) 



where p(A) is the spectral density of H = —ilb^ -\-V{x). This integral is convergent if 
a{\) ~ sgn (A) for large |A|; no cut-of is necessary. The classical static solution for M is 
then 

M=-eo + J ^(A)p(A)^. (52) 

We proposed the ansatz, 



where, 



777/ 

V{x) = mjo + -^-rile - 7o) (53) 



77(2;) = 1 + e-"", zy = 2msin^, ^ Q~\0)'yoQ{0), Q{e) ^ e'^'^/^ (54) 



Then, we showed that 

Q ^ 

(t(A) = sgn (A) for |A| > m, a{mcos-) — (55) 

2 71 

satisfies the self-consistency relation. (The value of cr(A) can be arbitrarily chosen for 
other values of A, since p(A) vanishes there.) Thus we have a one-parameter family of 
static solutions of the classical theory. The parameter ^ is a topological number that 
determines the behaviour of the solution as the center of mass variable goes to infinity. 
When ^ = 0, we recover the vacuum solution, V{x) = 77170. 

Although we have an explicit expression for p(A; x, y), it appears cumbersome to eval- 
uate the integral over A for M{x,y). What we mean by an exact solution is the above 
integral representation for M. 

Next we find the energy of this configuration. The integral for energy appears to be 
too hard to be calculated directly. We will instead use an indirect method. The variation of 
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E{M) under an infinitesimal change is tTH{M)6M. At a static solution, E{M) is invariant 
under all infinitesimal variations of M satisfying the constraint but that do not change its 
boundary conditions. However a change of $ in the above ansatz will change the behaviour 
of M at infinity, and E{M) need not be stationary with respect to this variation. ( 6 is 
a topological quantum number of the soliton.) We will get a simple expression for the 
derivative of energy with respect to 6. 

Let us calculate therefore 



Now we use the formula. 



/ 



d\ ,,,dp(X) f dXdX' a(X) - a(X') ^^^dH 



which follows from first order perturbation theory. The second term in the expression for 
^ becomes, 

f dXdX' a{X) - a{X') ,,,,,, ,^,dH 
' = J A-A- ' MA-)gMA)^. (87) 

The trace will be non-zero only if A = A', since H commutes with p(A'). Then we will 
have 

fdXdaiX) „ ^^^dH 

Now (7(A) is a constant on the continuum, so there is no contribution from it to this integral. 
As for the bound state, recall that ( for 6* 7^ 0, 2n) it is always separated by a gap from 
the continuum. We can continue <7(A) into this gap arbitrarily since p(A) is zero there; we 
could for example choose a{X) to be constant in some interval containing A^. Thus we see 
that / = 0. This can probably also be seen by more direct but tedious calculations. We 
have, 

dE{e) _ fdXda{X) _1 



dO 



Thus we find ( recall that the configuration with ^ = is the vacuum which has zero 
energy) 

£;(^) = -msin^. (60) 

TT 2 



16 



This agrees with the large N hmit of the spectrum of the Thirring model as obtained by 
the Bethe ansatz method, if we identify fx = ^m. 

By similar arguments we can also obtain the derivative of the energy density, 

dE(x,e) 2 6* ,t, , , , , 

— 00 — = -"^cos -iphix)i^B{x) (61) 

which leads to 

11^ 9 \ 9 
E(x,9) — — [— sin - tanhfmxsin -) r- logfcoshfmx sin -))] (62) 

This energy density is peaked around the origin and vanishes exponentially at infinity. 
This is precisely what we expect for a soliton. It should also be clear that the soliton 
quantum number 9 corresponds to the (large limit of) the discrete conserved quantum 
number of the fermionic theory. ( However, the 'baryon number' B = — trM is zero for all 
these configurations. This can be verified by using the regularization methods of [20]. 
Thus the situation is the opposite of that in 2dqcd; it is the abelian part of the symmetry 
that is 'confined'.) 

There are many additional properties of these solitons that can be studied. For ex- 
ample it would be interesting to obtain the time dependent solutions that represent the 
classical scattering of several of them. It would also be interesting to generalize the inverse 
scattering methods of Ref. [19] to infinite component classical field theories such as ours. 
Furthermore, it should be possible to generalize our results to Thirring models with several 
flavors. 

We thank T. Turgut for useful discussions. This work was supported in part by by 
the US Department of Energy, Grant No. DE-FG02-91ER40685. 
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